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Abstract Magnetic field behaviour in a spherically-symmetric accretion flow for parameters typical of single black holes in 
the Galaxy is discussed. It is shown that in the majority of Galaxy volume, accretion onto single stellar-mass black holes will be 
spherical and have a low accretion rate (10~ 6 - 10~ 9 of the Eddington rate). An analysis of plasma internal energy growth during 
the infall is performed. Adiabatic heating of collisionless accretion flow due to magnetic adiabatic invariant conservation is 
25% more efficient than in the standard non-magnetized gas case. It is shown that magnetic field line reconnections in discrete 
current sheets lead to significant nonthermal electron component formation. In a framework of quasi-diffusion acceleration, 
the "energy-radius" electron distribution is computed and the function describing the shape of synchrotron radiation spectrum 
is constructed. It is shown that nonthermal electron emission leads to formation of a hard (UV, X-ray, up to gamma), highly 
variable spectral component in addition to the standard synchrotron optical component first derived by Shvartsman generated 
by thermal electrons in the magnetic field of accretion flow. For typical interstellar medium parameters, a black hole at 100 
pc distance will be a 1 6-25 m optical source coinciding with the highly variable bright X-ray counterpart, while the variable 
component of optical emission will be about 18-27 m . The typical time scale of the variability is 10~ 4 sec, with relative flare 
amplitudes of 0.2-6% in various spectral bands. Possible applications of these results to the problem of search for single black 
holes are discussed. 
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1. Introduction 

Even though more than 60 years have passed since the the- 
oretical prediction of black holes as an astrophysical objects 
(Oppenheimer & Snyder 1939 1 in some sense they have not 
been discovered yet. To identify an object as a black hole, one 
needs to show that its mass exceeds 3M , its size is close to 
r g - IGM/c 1 and it has an event horizon instead of a normal 
surface - the distinguishing property of black holes which sep- 
arates them from massive compact objects of finite size in some 
theories of gravity (Will ll99'8l . However, only the two former 
criteria are used now for selection of black hole candidates of 
two types: a) with masses of 5-18 M Q , in X-ray binaries (see, 
for example, Greiner et al. 2001 1; and b) supermassive black 
holes in galaxy nuclei with masses of 10 6 - 1O'°M (Shields 
1999 1. Existence of the event horizon in such objects is usu- 
ally implied by the absence of periodic pulsations of the X-ray 
emission from strong regular magnetic fields (the black hole 
"no-hair" theorem) and I type X-ray flares due to thermonu- 
clear bursts of the accreted matter on the surface of the neu- 
tron star. At the same time, typical masses of X-ray pulsars 
and bursters are close to the typical neutron star value of 1.4 
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M while black hole candidates, missing pulsations and X-ray 
flares, have masses of 5-18 M (Miller et al. 119981 . The ab- 
sence of an event horizon in low-mass objects is not a proof of 
its existence in higher-mass ones. 

High accretion rates in X-ray binaries and active galactic 
nuclei result in the screening of regions close to the event hori- 
zon, and the most luminous parts of accretion flow are situated 
at distances of 10 - 100r ? (Chakrabarti 1996; Cherepashchuk 
2003 ) where general relativity effects are negligible. 

There is a very effective way to get information about the 
innermost parts of accretion disks in X-ray binaries as well as 
AGNs - the investigation of the broad (and sharp) iron K a fluo- 
rescent emission line (see review by Reynolds & Nowak 2003 1. 
Its intensity and shape depend on the accreted plasma distribu- 
tion and behaviour until the last stable orbit (0.62 r g for a ex- 
tremely spinning Kerr black hole and 3 r g for a Schwarzschild 
one) (Miller et al. 120041 Miniutti et al. 2004 1. However since 
the photons generated at the different distances from the hori- 
zon are mixed in the line profile it is not possible to extract the 
manifestations of gravitational fields close to the horizon only. 
This may be possible by study of variability of the iron line 
(Reynolds et al. l2"0Tj4l 
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At the same time, single stellar-mass black holes, which ac- 
crete interstellar medium of low density (10~ 2 - 1cm" 3 ), are 
the ideal case for detection and study of the event horizon. 
Shvartsman ( 1971 1 first demonstrated that an emitting halo of 
accreted matter forms around such objects and generates opti- 
cal featureless emission. The majority of such emission comes 
from the regions near the horizon at (3 5)r g . Spherical ac- 
cretion onto the single stellar-mass black holes has been studied 
in detail in the works of several authors (Bisnovatyi-Kogan & 
Ruzmaikin fT974l Meszaros [T9751 Ipser & Price [19771 1T9821 
and the main conclusions of Shvartsman have been confirmed. 

The most striking property of the accretion flow onto the 
single black hole is its inhomogeneity - the clots of plasma 
act as a probe testing the space-time properties near the hori- 
zon. The characteristic timescale of emission variability is r v ~ 
r g /c ~ 1CT 4 - 1CT 5 sec and such short stochastic variability 
may be considered as a distinctive property of black hole as 
the smallest possible physical object with a given mass. Its pa- 
rameters - spectra, energy distribution and light curves - carry 
important information on space-time properties of the horizon 
(Beskin & Shvartsman 1976 1. 

The general observational appearance of a single stellar- 
mass black hole at typical interstellar medium densities is the 
same as other optical objects without spectral lines - DC- 
dwarfs and ROCOSes (Radio Objects with Continuous Optical 
Spectra, a subclass of blazars) (Beskin & Mitronova 119911 
Pustilnik [T9T71 Shvartsman [T9T71 Beskin et al. EUUOl . The sug- 
gestion that isolated BHs can be among them is the basis of 
the observational programme of search for isolated stellar-mass 
black holes - MANIA (Multichannel Analysis of Nanosecond 
Intensity Alterations). It uses photometric observations of can- 
didate objects with high time resolution, special hardware and 
data analysis methods (Shvartsman 1977, Beskin et al. 119971 . 

In observations using the 6-meter telescope of the Special 
Astrophysical Observatory of 40 DC-dwarfs and ROCOSes, 
only upper limits for variability levels of 20% - 5% on the 
timescales of 10" 6 - 10 sec, respectively, were obtained, i.e. 
BHs were not detected (Shvartsman et al. U989al[T989bl Beskin 
et al.EDOOi. 

Recently, some evidences appeared that single stellar-mass 
black holes may be found among the stationary unidentified 
gamma-ray sources (Gehrels et al. 2002), gravitational lenses 
causing long-lasting MACHO events (Bennett et al. 2001) and 
white dwarf - black hole binaries detected by means of self- 
microlensing flashes (Beskin & Tuntsov 2002 1. In the last case, 
the mass transfer from the white dwarf is absent and a black 
hole behaves as a single one. 

In the present work we study spherical accretion onto a sin- 
gle stellar-mass black hole at low accretion rates of 10 8 - 10 13 
g/s. This corresponds to the range of interstellar medium den- 
sities of 0.002 - 0.1cm" 3 and a 10 M Q object moving with a 
velocity of 20-40 km/s (Bondi & Hovle ll944l . This is true for 
about 90% of the Galaxy volume (McKee & Ostriker fT977l l. 

Spherical accretion with equipartition of energies, i.e. with 
roughly equal densities of the magnetic and kinetic energies of 
plasma has been considered in many papers ( Shvartsman l 1 97 1 1 
Bisnovatyi-Kogan & Ruzmaikin 119741 Kowalenko & Melia 
119991 Ipser & Price 119771 IT982). The uniqueness of our ap- 



proach is in taking into account the significantly non-thermal 
nature of electron energy distribution function (its synchrotron 
emission determines the appearance of a black hole). It may be 
roughly considered as a superposition of two components (this 
approach is known as "hybrid plasma", see Coppi ( 1999) and 
references therein) - thermal electrons and accelerated electron 
beams, formed in current sheets where magnetic energy is dis- 
sipated in a way similar to solar flares (Pustilnik ll978l [T997 1. 
The latter process supports the equipartition of energies. As a 
result, the emission of the accretion flow consists of a quasi- 
stationary "thermal" part with a wide-band spectrum from in- 
frared to ultraviolet, and a highly variable flaring nonthermal 
component. Each such flare is generated due to the motion of 
the accelerated electron beam in the magnetic field. Its light 
curve carries information on the magnetic and gravitational 
field structure near the black hole horizon. Nonthermal lumi- 
nosity reaches several percents of the total luminosity and may 
even exceeds it at low rates, while its spectrum covers spectral 
bands from optical to hard X-ray. This result leads to possible 
modifications of the search strategy. 

In §2, the main characteristics of accretion onto single 
black holes in the Galaxy are discussed. 

In §3, the electron distribution function in phase space is 
built, in §4, the thermal and nonthermal component luminosi- 
ties are determined, and in §5 the shape of its spectra is studied. 

In §6, the temporal behaviour of single electron beam emis- 
sion is studied and some conclusions on the variability of ac- 
cretion flow are made. 

In §7, the main results of this work are summarized, and in 
§9, possible directions of future work are discussed. 

2. Nature of accretion and basic parameters of the 
model 

2.1. Accretion rate for different parameters of black 
holes and interstellar medium 

Contemporary models of massive star evolution, its dynam- 
ics and data on black hole candidates in X-ray binaries and 
microlensing events suggest that the most probable mass M 
of a single black hole is 10M o (Greiner et al. 2001 Fryer & 
Kalogera 2001 Agol & Kamionkowski 2002 1, while its veloc- 
ity is in the 10-50 km/s range and the gas capture cross-section 
is defined by the Bondi radius r c (Bondi & Hoyle 1944 1: 



where c s is the sound speed in the interstellar medium. 

The interstellar medium consists of at least three compo- 
nents (McKee & Ostriker fT977t - cold (T ~ 10 2 K) and dense 
(n ~ 10 2 cm" 3 ) neutral hydrogen clouds, warm partly ionized 
hydrogen (n ~ 0.1 - 0.2 cm" 3 , T ~ 10 4 K) ones and fully ion- 
ized coronal hydrogen (n ~ 0.002cm" 3 , T ~ 10 6 K) clouds. 
Motion may be subsonic as well as supersonic - the sound 
speed changes in the 1.5 — 150 km/s range. 

For the black hole moving supersonically in a uniform 
medium, an accretion rate is determined by the behaviour of 
matter behind the gravitation centre, where the tangential gas 
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Figure 1. Dimensionless accretion rates m = Mc 2 IL e dd as a 
function of interstellar medium density for various black hole 
velocities. 

velocity component vanishes in the shock wave and the gas 
falls towards the black hole in a wide cone(Bondi & Hoyle 
[T944I Shvartsman fT97Tl Illarionov & Sunvaev [19751 Font & 
Ibanez 1998 1. For the typical case of interstellar medium that 
is collisionless on scales of the capture radius, a mode of ac- 
cretion from a very thin "tail", where the tangential velocity 
vanishes due to the magnetic field, is possible. The thickness 
of such a tail may be estimated as 



rug ( B 2 J8tt \ 
r c ~\pV 2 /2) 



1/2 



V" 



(2) 



where B^ is the magnetic field strength at infinity, V is 
the black hole velocity and rough equipartition of thermal and 
magnetic energies in the interstellar medium is taken into ac- 
count. Accretion flow is continuous due to interstellar magnetic 

fields - the proton Larmor radius is 10 7 10 8 cm which is at 

least 10 6 times smaller than r c . The interstellar magnetic field 
in the Galaxy (10 4 — 10" 5 G) is frozen-in. The accretion rate 
for such a regime is determined by the expression (Bondi & 
HoyleEHUl 



M = 



4nG 2 M 2 p 

( y 2 + c 2 )3 /2 



(3) 



where p in the interstellar medium density. 

Using the expression for the Eddington accretion rate 
Medd = L edd lc 2 = = 1.4 ■ 10 18 (M/10M o ) g/s, where m p 

is the proton mass, cr T is the Tompson cross-section, we have 
a useful normalization, measuring the mass in units of 10M o , 
density in 1 cirT 3 , and velocity in 16 km/s 



m = MIM edd = 1.3 ■ 10- 5 M Ioni (V 2 + c 2 s )7 3 ' 2 



(4) 



In Fig [2 the dependence of m on n is shown for various ve- 
locities of the black hole motion (the dependence of the sound 
speed c s on the temperature T due to hydrostatic equilibrium 
of the ISM is also taken into account). Only in the cold clouds 
the accretion rate may reach the Eddington level, but due to 



the relative rareness of such clouds (5% of the Galactic vol- 
ume) this case is improbable. Note however that a black hole 
may initially be born in such a cloud, and so become a very 

bright source with a luminosity up to 10 38 10 40 erg/s (and 

so, be a so-called "ultraluminous" source (Roberts et al. 2002 1, 
like those observed in other galaxies). For the warm hydro- 
gen, whose volume fraction ~ 40%, the accretion rate will be 
10~ 6 - 10~ 8 , while for hot hydrogen (~ 50% volume) it will 
be significantly lower. Later, we will assume such cases as the 
most typical for single black holes in the Galaxy. 

2.2. The role of interstellar medium inhomogeneities. 
Spherical or disc-like accretion? 

The observational appearance of the accreting black holes de- 
pends crucially on the regime of matter flow near the relativistic 
object. So, in binaries the captured angular momentum of mat- 
ter is high enough, and a disc-like accretion regime takes place 
(Shakura & Sunvaev ll973> . For the case of accretion from in- 
terstellar medium, however, this problem is more complicated. 

In the recent series of papers on the analysis of accretion 
onto single black holes, it was reasoned that a disc-like regime 
is realized for such a case also (Fujita et al. [T9981 Agol & 
Kamionowski 2002 Chisholm et al 2003 1. However, it seems 
they used overestimated values for the captured angular mo- 
mentum. 

Davies & Pringle (1980 1 consider analytically the prob- 
lem of accretion onto the moving gravitational centre of non- 
uniform interstellar medium (inhomogeneities of density as 
well as of velocity are considered) and show that, in the first 
order of inhomogeneities, the magnitude accretion rate is still 
described by the Bondi-Hoyle formula Q, and the captured 
angular momentum is zero. 

The latter result may easily be understood by considering 
a simple model of accretion from the tail in which the traver- 
sal velocity component vanishes. This tail, on which the angu- 
lar momentum is zero by definition, will not, in general, be a 
straight line, and the capture cross-section will not be a circle, 
and so, the total captured angular momentum will become zero 
too. In a more realistic picture, however, the role of gas pres- 
sure is important, the matter accretes in a wide cone, and the 
total captured angular momentum may be significant, although, 
as it was shown by numerical simulation results of Sawada et 
al. (1989), Ruffert STW1\ [T999>. it is much smaller than the 



(5) 



usually-used estimation of Illarionov & Sunyaev ( 1975 1: 
L = \BV K r c , 



(where B = — , 4r — relative density and velocity variations 



on 2r c scale, Vk 



V 2 +c? 



Keplerian speed at the capture 



radius). So, this estimation may be considered as an upper limit 
in a very narrow class of objects. 

Detailed analysis and numerical simulations leads to the 
more realistic expression for the captured angular momentum 
(Ruffert fT997l fT999i: 



Z~0.1L 



(6) 
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For the spherical accretion regime the specific angular momen- 
tum of the captured matter must be smaller than that of the 
black hole last stable orbit, i.e. / < y/3cr g for Schwarzschild 
metric. Then 



where a 2 as 1/3 (which corresponds to equal amounts of grav- 
itational energy transition into kinetic, magnetic and gravita- 
tional ones). Therefore, for the parameters of the accretion flow 
we have 



P < 40 Vr5 



(7) 



acR 



-1/2 



(12) 



Due to the characteristic dispersion scale in the interstellar 



medium being of the order of (A V) 2 = l.l(r/lpc) u - /& (^) 
(LarsonHMU Falgarone & Phillips fT990l . then 



krru2 P ~ 



Ipc s 
where Vq = yjV 2 + c\. Finally 



Vq > 17M? 138 km/s 



(8) 



(9) 



Therefore for nearly any black hole velocity, interstellar 
medium turbulent motion cannot prevent realization of the 
spherical accretion regime. 1 Furthermore, density fluctuations 
cannot prevent it either. For -j- ~ 
1995 1, it may be easily shown that 



(l^) 6 (Armstrong et al. 



Vo > 3.1 M^f km/s 



(10) 



It is clear that even in cold clouds, density fluctuations can- 
not lead to the disc accretion regime. 

2.3. Radial structure of the flow 

The general solution of hydrodynamical problem of accretion 
onto a non-moving gravitating center derived by Bondi ( 1952 1 
determines radial profiles of various accretion flow parameters 
as a function of the distance to the black hole, accretion rate 
and gas adiabatic index. It is important for us that for y < 5/3 
(that is always true for interstellar gas) there is a "sonic point" 
in the flow, passing which the gas motion becomes supersonic, 
and the gas velocity near the gravitating center has an asymp- 
totic behaviour v oc r~ 1/,z (and sop oc r~ 3 ^ 2 for density). In an 
approximate description, we may extrapolate such a behaviour 
to the capture radius scale and use it for the whole flow. 

In the black hole tail, where the matter flow stops, the 
thermal, gravitational and magnetic energy densities become 
nearly equal (since in the collisionless case, the matter is 
stopped due to the magnetic pressure and plasma oscilla- 
tions (Illarionov & Sunvaev tl975» . This equipartition of en- 
ergies (at least magnetic and gravitational ones) is preserved in 
the following infall (so-called Shvartsman equipartition theo- 
rem (Shvartsman 1971 1), so we may assume that (Bisnovatyi- 
Kogan & Ruzmaikin 1 1 974> 



B 1 



1 



8^ = 2 PV 



a 



,GMp 



(ID 



1 As already noticed, an obvious exception is the case of cold clouds 
where sonic speed is low and at v < 16 - 17 km/s a disk may form. 
We are planning to discuss this case in a separate paper. 



M 



Anr 2 v la-jar , 



■R 



-3/2 



B 1 i am„(rm - n 

— = - pv 2 = p - R~ 5/z , 

8n 2 \ajr g 



(13) 



(14) 



where the dimensionless values for radius R = r/r g and the 
accretion rate from equation (0J are used. Numerical values of 
these parameters are 



n = — = 4.33 ■ 10 12 m-sMTiR 
m 



lp-3/2 cm -3 



B = 8 • 10 4 th^Aq^R 



1 / 2 M- 1 /2i ? -5/4 



Gauss. 



(15) 
(16) 



The magnetic field has a quasi-radial sectorial structure (the 
radial component grows much faster than the tangential one, 
which is proportional to the square root of distance; magnetic 
field lines are stretched). 

The assumption introduced earlier on the equipartition of 
magnetic, kinetic and gravitational energy requires the exis- 
tence of magnetic flux (and, therefore magnetic energy) dissi- 
pation mechanisms. It was first noted by Shvartsman J197H : 
their possible observational appearances were discussed by 
Illarionov and Sunyaev dl975l l: a model-independent estima- 
tion of the energy dissipation rate was proposed by Bisnovatyi- 



Kogan & Ruzmaikin ( 1974 1 and from a different point of view 
by Meszaros ( 1975); alternative approaches were discussed by 
Scharlemann dl983t and Kowalenko & Melia J1999I . 

The dissipation rate of such a mechanism may be estimated 
as follows (Bisnovatyi-Kogan & Ruzmaikin I1974i . Scaling 
laws for the magnetic field in a given volume element for the 
frozen-in field and equipartition are correspondingly 



d B 2 \ 

dt 8n l frozen-in 
dt $7T J e q U ip ar titi on 



vB z 
5 vB 2 



B 2 



1 



8tt 2 PV 



,GMp 



(17) 



(18) 



(19) 



For the preservation of equipartition state, the power equal to 
the difference of these expressions must be dissipated in a vol- 
ume element: 



dE 



3 vB 2 



dVdt 2rSn 

or the same for a spherical shell 

dE 3a 2 Mc 2 



dRdt 



4 R 2 ' 



(20) 



(21) 
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which, by integrating over the whole volume, gives 



dE 



1 



— = -a z Mc z = -Mc\ 
dt 4 4 



(22) 



which means that in the case of equipartition for the spheri- 
cal flow described earlier, as much as 25% of infalling matter 
rest energy is released through this dissipation mechanism only. 
This leads to an additional super-adiabatic heating of the gas 
(Bisnovatyi-Kogan & Ruzmaikin ll974l Meszaros ll975l Ipser 
& Price l 19771 . and unavoidably changes the temperature radial 
profile and the luminosity of the accretion flow. 

2.4. Flaring dissipation of the magnetic field and 
electron acceleration in the current sheets 

It is clear from simple physical reasons that dissipation of the 
magnetic field means that it is no longer frozen-in, i.e. the mean 
conductivity becomes much lower and relative motion of the 
magnetic field and plasma appears (the mean diffusion time be- 
comes comparable to the free-fall one), currents begin to flow 
and heat the gas. These processes may take place either con- 
tinuously in turbulent accretion flow as an Ohmic dissipation 
of the magnetic field (Bisnovatyi-Kogan & Ruzmaikin 1 1 9741 
Bisnovatyi-Kogan & Lovelace 1997 1 or in compact enough 
separate regions. In the latter case, the dissipation process has a 
highly discrete nature (this possibility was noted by Bisnovatyi- 
Kogan & Lovelace (2000l) of single events of magnetic field 
line reconnections. 

In the present work we consider the latter case. It seems 
to be more realistic due to the fact that Alfven velocity which 
determines the speed of energy exchange between such a re- 
connection region and the surrounding plasma is nearly equal 
to the free-fall velocity at which a high magnetic field gradient 
forms. This is analogous to the case of continuous energy sup- 
ply to the magnetic field inhomogeneities in the solar corona 
and its flaring (discrete) dissipation and conversion to acceler- 
ation of particles and anomalous turbulent heating. This is the 
case of formation of turbulent current sheets, which leads to 
fast magnetic field line reconnection, avalanche-like growth of 
energy release (the flare itself) and threshold switch-off of en- 



ergy dissipation processes (Sweet 1969; Petchek 1964; Parker 



1979J Spitzer fT954l SyrovatskiiQ981 ). Multi-frequency obser- 
vations of solar flares support this mechanism. On the other 
hand, the similarity of statistical properties of flaring activity of 
the Sun (Lu & Hamilton [T99TI Lu et al. [T993l and UV Cet stars 
(Gershberg 1989 1, X-ray binaries (Kawaguchi & Mineshige 
11999 1 gamma-ray bursts and active galactic nuclei argues in 
favor of the universality of such processes. So, we may assume 
that the main mechanism providing magnetic energy dissipa- 
tion in the accretion flow is a reconnection of the magnetic field 
lines in the current sheets (regions of high magnetic field gra- 
dient). 

Energy dissipation in the current sheet itself may be con- 
sidered as a simple Joule heating Q = jE — j 2 /cr, where 
j = (c/4-n) ■ rotB oc — - current density, a - current sheet thick- 
ness, cr - some effective conductivity. When current reaches 
some threshold value, j c = en c Cj, where n c - electron num- 
ber density, c, - ion sound speed, it generates strong ion- 



acoustic turbulence that lowers the conductivity by 9-10 orders 
of magnitude (electrons begin to collide with plasmons and 
heat them). The analysis of particle acceleration in such a cur- 
rent sheet (which may be considered as a superposition of di- 
rect electric field acceleration and diffusion - elastic scattering 
on plasmons) was performed in Pustilnik ( 19781 1199*71 . In the 
framework of this model, the accelerated particle energy distri- 
bution (electrons mostly) naturally appears to be a power-law 
up to energies large enough (Pustilnik 1978 1 and has a shape 



Mr) 



i/r 



(23) 



where F corresponds to the mean energy of a particle in a beam. 

We estimate the maximal energy a electron reaches by ac- 
celeration in such current sheet. By using expressions for gas 
density and magnetic field strength J13i . <ll4i for low accretion 
rates (m < 1CT 5 ), we have for the event horizon n < 4 • 10 12 
cm -3 and B g < 1.7 • 10 5 Gs. Near r g velocities of electrons 
and plasmons are nearly equal to the speed of light, and so, the 
maximum current density is j max = enc, and the maximum ef- 
fective field is E max - p?, where cr* is anomalous resistivity; 
for our case cr* ~ ^tj oe , where u oe - Lengmure frequency 
(cj oe = 5.65 • 10 4 • n 1 / 2 ). Therefore, the maximum gamma-factor 
is 



e 2 nr g 



cr*m e c 4.5 • I0 5 m e c 



10 s 



(24) 



It is very difficult to estimate the fraction £ of the dissi- 
pated energy carried away by the accelerated particles - elec- 
trons during the acceleration process generate ion-acoustic and 
Lengmure plasma oscillations. Also, topology changes in re- 
connections lead to global matter motions. A study of such pro- 
cesses in the solar flares shows that particles (electrons mostly) 
carry away from 10% to 50% of energy stored in the mag- 
netic field inhomogeneities (Hudson & Rvan ll99"5t . We assume 
£ = 0. 1 as a reasonable lower limit for further estimations. The 
(1 - £) fraction of the dissipated energy goes into the surround- 
ing plasma heating. In the collisionless case this is due to gen- 
eration of MHD turbulence in the current sheet. Also, as the 
speed of plasmons is nearly equal to the Alfven speed, which 
in turn is equal to the free-fall velocity for the equipartition 
case, the heating is nearly uniform through the whole accretion 
flow. It is also difficult to determine the motion of the acceler- 
ated particle beam. Of course, electrons are moving along the 
magnetic field lines, tracing its topology. Its motion is local due 
to the relative smallness of Larmor radius in comparison to the 
characteristic scale of the accretion flow, Schwarzschild radius 



y 2 m e 



2oj_ 
hurl 



« 1 



(25) 



For a typical single stellar-mass black hole accretion rates, 
electron-electron and electron-ion energy exchange is highly 
ineffective (this question has been discussed in detail in 
the framework of advective disk models, see Mahadevan & 
Quataert ( 1997 1 and Bisnovatyi-Kogan & Lovelace (2001 1), so 
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trum of flares on various objects - from active galactic nuclei 
to UV-Cet stars in the framework of a single mechanism. 

To complete the picture, we briefly review other models. 
Fast reconnections leading to particle acceleration were con- 
sidered by Lazarian & Vishniac ( 1999 2000), the role of var- 
ious instabilities in creation of anomalous resistivity and par- 
ticles acceleration has been considered by Birk et al. {Tjj99 1, 
and processes of particle acceleration in various astrophysical 
objects by Bisnovatyi-Kogan & Lovelace ( fT997l l2000l 1200 1 i . 

The dissipated magnetic energy is converted mostly to ac- 
celeration of electrons. In the framework of our paper this 
means that for analysis of emission processes we may take into 
account the electron component only. 



Figure 2. A schematic picture of the accretion flow structure. 
Reconnection regions and ejected plasma flows are shown. The 
direction of the particles acceleration (current flow) is perpen- 
dicular to the viewing plane. 



we may neglect collisional energy losses for the accelerated 
particles. 

Furthermore, in the case of a Maxwellian distribution of 
the background plasma electrons and moderate mean gamma- 
factors of the accelerated nonthermal electrons (10 - 10 2 of 
the mean thermal one), the total electron distribution f(y) al- 
ways (for energies above thermal peak) satisfies the inequal- 
ity < 0, and so is stable to the generation of plasmons 
with wave vectors parallel to the motion direction (Kaplan & 
Tsvtovich ll973> . For the same reason, all generated instabil- 
ities vanish rapidly and the beam is stabilized. However, the 
induced generation of plasmons with non-collinear wave vec- 
tors is possible; this leads to the dissipation of the beam with a 
characteristic time of tdis ~ lO 2 -^^^, where n* is the electron 
density of the beam (Kaplan, Tsvtovich ll973> . It is easily seen 
that 7^ ~ r~ 3 ' 4 , so, the beam born near r„ has more chance of 

surviving. Also it is seen that for jf- > 10 2 the beam has not 
enough time to dissipate in the free-fall time scale, so, we may 
neglect this scattering and thermalization of the beam and con- 
sider it as a system of non-interacting electrons in the external 
magnetic field. 

The magnetic field dissipation in the current sheet is ac- 
companied by significant dynamical effects; the reconnected 
magnetic field lines are ejected through the current sheet ends 
carrying away the gas that lowers the current sheet density (see 
Fig-E}. 

The results of the solar corona and chromosphere observa- 
tions (Dere 1996 Innes et al. H997i show that the reconnection 
speed is often significantly less than the Alfven one (« 0.1 V^). 
In such a case the density and temperature fluctuations may not 
have enough time to vanish, and this may lead to additional 
variability of the accretion flow. 

The magnetic energy dissipation mechanism proposed here 
is not a unique one. However, it is in a good agreement with the 
observations of the Sun (Pustilnik ll997> . Also, it provides an 
opportunity for interpretation of the universal energetic spec- 



3. Derivation of electron distribution 

3.1. Note on adiabatic heating 

At the low accretion rates considered, the gas is fully colli- 
sionless, and so the particle energy distribution is determined 
by superposition of three different factors - adiabatic heating, 
synchrotron cooling and nonthermal heating due to particle ac- 
celeration in the current sheets. The former dominates at large 
distances, and so its precise treatment is very important. 

At the level of single particle motion, adiabatic heating of 
collisionless magnetized gas is due to the conservation of the 
adiabatic invariant (Landau & Lifshitz 1971 1 



/ = 



2eB 



(26) 



of the charged particle in the magnetic field evolving slowly on 
a timescale of a single Larmor revolution. This basically means 
the conservation of the phase volume element per particle p] ■ 
V = const. Here we neglect the motion parallel to the magnetic 
field line - only the perpendicular momentum p, grows due to 
compression, while the parallel one p\\ is determined mostly by 
the initial conditions. 

Note that classical relation between the spatial part of the 
phase volume and the gas density Vp - const, that is true for 
normal "gas in a box" of classical thermodynamics, is gener- 
ally wrong for magnetized plasma. An exception is the case of 
isotropic gas compression with B oc l~ 2 oc p 1 ^ - in this case the 
magnetic field acts as the walls of a box, and this leads to the 
usual equations of state e oc p 5 ^ 3 for the non-relativistic case 
or e oc p 4 / 3 for relativistic one. Here we use the proportional- 
ity of the particle energy to its perpendicular momentum, that 
is true for two opposite cases - either for p\\ oc p, (effective 
isotropization) or for p\\ <k p, (collisionless case, only perpen- 
dicular momentum increases). 

In the case of the accretion flow, however, the gas com- 
pression is significantly anisotropic (the matter element is even 
stretched in the radial direction proportional to r ^ 2 ), and the 
magnetic field itself is not perfectly frozen-in either. So, the re- 
lation of the phase volume element spatial part to the density 
becomes 



dV 3^' 3 dZ? 15 dr 5 dp 5 dn 
~V~~ ~p7 ~ ~2~B ~ _ ~8~T ~ ~4~p~ ~ ~4~ ' 



(27) 
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where the radial dependencies (I13M14> of the accretion flow 
parameters are used. 

Note that in contrary to assumptions made in most articles 
on this subject (Bisnovatyi-Kogan & Ruzmaikin l 19741 Shapiro 
I1973bl Meszaros lT975l Ipser & Price 119771119521 Mahadevan 
& Quataert 1997 1 we cannot use the first law of thermodynam- 
ics in the form of "energy per particle" 



\nl \n 



(28) 



(the dots here represent the contribution of non-adiabatic pro- 
cesses) to describe the relation of the particle mean energy 
and the gas density as the density changes don't reflect the be- 
haviour of the particle "walls". A correct form of this equation 
must take into account the fact that the number of particles N 
inside the spatial part of each particle conserving phase volume 
element is no longer constant 



d(eV) = -pdV + -dN + 



(29) 



Using 1271 . the change of the number of particles may be ex- 
pressed as 



dV d(l/n)\ 1 

dN = nV\ — - = — Vdn. 

\ V (In) ) 4 



(30) 



For relativistic gas the equation of state has the form p = e/3, 
and so the single particle energy evolution is described by 



(31) 



where dots represent the possible contribution of non-adiabatic 
processes per particle (thus, such particles behave like a gas 
with specific heat ratio 11/6). 

For non-relativistic gas in a similar manner (p = §e) 

d $=-H$ + -' (32) 

like the gas with a specific heat ratio of 17/12. 

Equations (1321 and ( 13 1 1 may be rewritten in a form useful 
for comparison with the incorrect expression ( 1281 : 



d (^) = - 5 4 Pd {l 



+ ... . 



(33) 



It is easily seen that correct consideration of adiabatic heat- 
ing makes it 25% more effective than in the case of ideal non- 



magnetized gas accretion of Bondi ( 1952 1, that has a large in- 
fluence on luminosity and spectral shape of the accretion flow 
emission. 

3.2. Radial temperature distribution 

Note that adiabatic heating alone does not change the shape of 
the particle momentum distribution, so initially thermal distri- 
bution always stays thermal. 

Current sheet spatial scales are usually much smaller than 
the whole accretion flow one and, therefore, the fraction of 
the accelerated particles is small, and so the total (significantly 



nonthermal) electron distribution may be considered as a su- 
perposition of the purely thermal one for the background flow 
particles and purely nonthermal for the ones accelerated in the 
current sheets (this is the approach known as "hybrid plasma" 
of Coppi ( 1999 1). Also we assume that for low accretion rates 
the non-adiabatic heating and radiative energy losses do not 
change the shape of the thermal component distribution (while 
changing its mean energy) 



f(R,y) = f t (R,y) + H nt (R,y). 



(34) 



Note that this distribution is not normalized, and only its shape 
has physical meaning. So, for example, the ratio of nonthermal 
to thermal electron densities at some radius R may be expressed 



Hg(g) = (fnt(R) 

n,(R) MR) 



(35) 



where f„ t (R) and /,(/?) are integrals of the corresponding distri- 
bution functions over the range of y. 

Thermal particle distribution function may be written as 



(36) 



where the usual dimensionless expression for temperature r = 
kT/m e c 2 is used. This gives a Maxwellian local energy distri- 
bution and radial density slope p <x R"^l 2 . 

The temperature distribution t(R) may be determined by 
solving the energy balance equation taking into account heat- 
ing due to adiabatic compression and magnetic field dissipation 
and radiative losses. Note that electrons become relativistic at 
some large radius R re i, while protons remain non-relativistic 
until the event horizon (Bisnovatyi-Kogan & Ruzmaikin 1974 1, 
and so their heating rates differ by a factor of 2 ( 13 1 132b . In the 
case of ineffective energy exchange (at low accretion rates col- 
lisions are very rare; other mechanisms noted in Mahadevan 
& Quataert ( 1997) are also ineffective); this could have led to 
electrons much colder than protons, but due to preferred heat- 
ing of electrons by non-adiabatic processes (Bisnovatyi-Kogan 
& Lovelace 2000 2001 1 its energies always remain roughly of 
the same order of magnitude. Of course their gamma-factors 
will differ by about 40 times, and so the main contribution to 
the accretion flow radiation is due to electrons. Therefore, only 
electron temperature is interesting for us. 

As it has been noticed before d20l >. the non-adiabatic heat- 
ing rate may be expressed as 



dE 3 v B 2 

* = (1 -0Sv* 88(1 -0275? ■ 



(37) 



The main mechanism of radiative losses at low accre- 
tion rates is synchrotron radiation. Its rate may be written as 
(Lightman & Rvbicki [T979l 



4 —B 2 

For a Maxwellian distribution 

12 1 -^y 2 = 12r 2 

m e c l 



(38) 



(39) 
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Taking into account non-adiabatic terms, the energy bal- 
ance equation per particle (I32> may be written for a non- 
relativistic region of the flow (R > R re i) as 



d e 5 £ 1 A-O 



dr n 



Anr v n 



(40) 



For a non-relativistic electron gas e = ^nkT; also we may 
neglect energy losses and rewrite it in variables R and t as 



dr 5 t 
dR ~ ~~4R 



a m 



P d-2 



2 m, 



(41) 



This gives the value of the radius for r = 1, R re i w 6400 for 
£ = 0.1, that significantly exceeds R re i as 1300 in the ideal 
gas approximation used by Bisnovatyi-Kogan & Ruzmaikin 

In deeper regions, electrons become relativistic, their en- 
ergy density is e = 3p = 3nkT, and so for R < R re i the energy 
balance equation (13 1 1 combined with I13i . J 1 4b and @ gives 



dr 
dR 



5 T 
"8R 



a 2 m„ . 4 m„ mr 2 
4 m„ 3 m P k z 



(42) 



The boundary condition is T(R re i) = 1. 

An analytical solution of this equation in general is very 
difficult, but for low accretion rates we may neglect the influ- 
ence of radiative losses and get a solution in the form 

2a 2 m„ . / 2a 2 m p \ I R re i\ 5,& 

t(R) = 1 -(1 -0 — '-)[-£) -(43) 

3 m e \ iRrei m e j \ R I 



This expression may be substituted into d36i to get the final 
expression for the thermal electron distribution. 

The applicability of this approximation may be estimated 
by comparing the timescales of the electron radiative energy 
losses and of the free-fall near the event horizon 



tff a m p 

— = -ym < 1 for m < 10 , 

trad 3 m e 



(44) 



where y x 100 is assumed. For higher accretion rates equation 
( 1421 may be easily solved numerically. 

3.3. Distribution function of the non-thermal 
component 

We can build an expression for nonthermal component distribu- 
tion function. We may assume that all electrons are relativistic 
as only these have significant observational appearances, and 
electrons are relativistic in the most interesting regions near 
the horizon (well inside the relativization radius R re i * 6000). 
So, the fraction of low-energy nonthermal particles is negligi- 
ble, and nearly all electrons have a "lifetime" not less than the 
characteristic free-fall time scale (see Sec. l2.4l i. 

So, we may neglect nonthermal particle interactions with 
electrons and plasmons of the background flow (see discussion 
in Sec. 12. 4t and assume that they evolve due to adiabatic heat- 
ing and synchrotron energy losses only. Note that this leads 
to the situation where thermal electron energy grows slightly 
faster than the nonthermal one due to additional heating by 
plasma oscillations ejected by current sheets. 



The energy evolution of single nonthermal electron is de- 
scribed by 



dy _ 1 m p y- 5 y 
dR ~ 3^T m F ~ 8R' 



(45) 



where the first term corresponds to synchrotron losses and the 
second one to adiabatic heating. 

For the initial energy yo at Rq it has the solution 



r 



yo 



C 1 (R,R )y Q + C 2 {R,Ro) 
where 



(46) 



CdR,R ) 
A 

C 2 (R,R ) 



1 - 



R 
R~o 



13/8\ 



39 m e 



5/8 



The nonthermal component at some radius R consists of 
non-interacting electron beams generated at all radii Rq > R. 
The evolution of the distribution function of each such beam 
(assuming that the ejection process is stationary and the initial 
beam distribution function is fh(Ro,7o)) i s as follows: 



MR,y) = MRo,7o) 



dyo 
dy 



f(Ro,7o) 



(CiOR,flo)yo + C 2 OR,flo)) 2 



(47) 



(48) 



C 2 (R,Ro) 

The initial beam distribution has a form (see i23\ ) 

f m ^ MRo)(nRo)\ 3 I mo) 

MRo, yo) = tt— exp 

i(Ro) \ yo / \ yo 

where the mean energy 



yo = r = At 

is assumed to be by a fixed factor A greater than the local ther- 
mal component one and fi,(R) describes the radial distribution 
of the accelerated particle ejection rate dN/dt. The latter may 
be computed using (12 1 1 as 



1 gdE 6na m p r g c fyh 



dN 

dRdt ~ m e c 2 T(R) dRdt ~ ~ 4 m e o- T R 2 T(RY 
which gives for the total ejection rate 



dN 6na 2 m p r g c 



dt 



4 m e ctt 



£m 



CO 

s 



dR 



R 2 r(R) ' 



(49) 



(50) 



and the radial distribution of particle ejection 



MRo) 



fl() 



r(R )R 2 



where ao 



• oo 

V i 



dR 
T(R )Rl 



(51) 
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By integrating ( 1471 over all Rq > R we may get final ex- 
pression for nonthermal electron distribution which consists of 
all beam electrons ejected at greater distances 



f n ,{R,y) 



/aod 



a dR T(R ) [Ci(R,R )y + C 2 {R,Rq)] 2 



C 2 (R,Ro) 



x exp - 



Jo 



(52) 



Now we determine the f coefficient of the total electron 
distribution J34t . By combining (I35> with trivial expressions 
for the energy density ratio 

e nt (R) = n nt (R) y,(R) 

e,(R) n t {R)y n ,{R) 

and for the energy density of thermal electrons with a 
Maxwellian value y,{R) = 3t(R) 



e, — 3m e c T(R)n, 



we may get 
e„ t (R)^ 



i = 



m e c 2 n,(R) 



oo 

V 

v i 



yfn,(R,y)dy 



(54) 



(55) 



The fraction of nonthermal particles is always small, so we 
may replace here the thermal particle number density n t (r) with 
the total one n(R). 

The nonthermal electron energy density at R may be com- 
puted by considering an elementary spherical shell R ~ R + 6R 
and integrating the nonthermal energy release in it during its 
whole backtrace free-fall history 



£m(R) 



1 



4nr 3 g R 2 dR 



5RodR (jdE_\ y b (R) 
dRodtj T(Ro)' 



f SRpdRp 



(56) 



where the latter multiplicative term corresponds to amplifica- 
tion of the mean gamma-factor of each accelerated electron 
beam ejected at Ro in free-fall until R. This quantity may be 
written using the results of AppendixlAlas 



n(R) 

URo) 



5/8 



AWq) 



R 



1 - 



R 
R~o 



13/8\ 



(57) 



This combined with the expression for the local nonthermal 
energy dissipation rate Oil and with the scaling for the given 
spherical shell thickness at free-fall SRq °c Rq 
final expression for the ^ coefficient 



-1/2 



leads to the 



3 m„ , r - 

4 m P 



s 



y,f nt {R,y)dy 



V 1 



J R n 



5/8 



AT(R Q ) 



R 



1 - 



R 
R~o 



13/8\ 



(58) 



A sample fraction of the nonthermal to total densities for 
different accretion rates is shown in Fig[3]as a function of dis- 
tance from the black hole. This fraction is small. Fig|4]shows 
the shapes of thermal and non-thermal distributions for two dif- 
ferent distances at low accretion rates. 




71 = 10 
71=10" 
71=10 
71=10 
71=10 



Figure 3. Nonthermal electron density fraction n„i(R) / n(R) for 
an M - IQM Q black hole, £ = 0.1 and the number of accretion 
rates. 

4. Emission spectrum 

The emission spectrum of a single electron for an observer 



at infinity has the shape (Lightman & Rybicki 1979 Shapiro 
I1973a> 



COS 0* 



L Y — In Jv'— —dcosO, 



I 



(1 -ygCOSfl) 2 



(59) 



where relativistic effects of time contraction, gravitational red- 
shift, Doppler effect and the capture of some emission fraction 
by the event horizon are taken into account. The event horizon 
angular size for a free-falling emitter is 



27 1 

| COS 6»*| = J I r(l - -), 

4R 2 R 
where cos 9* < for R < 1.5. The quantity 

dr 1 v/c 



dtl-rJr I 2 , 2 



(y 2 /c 2 + 1 - r g /r) 



1/2 



represents the falling velocity of matter in the distant observer 
frame, and the frequency shift is given by 

, 1 - (v/c) cos 9 



V(l - v 2 /c 2 )(l - l/R) 



By substituting synchrotron emissivity from Lightman & 
Rybicki (1979) we get 



V3e 3 fisin(/r /v' 
' 4-nm e c 2 \ v c 



F(x) = x J K 5/3 (&dt, 

X 

3y 2 eB simfr 



Vc = 



Anm e c. 



(60) 



(61) 



(62) 
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Thermal - 
Non-thermal 




Figure 4. Energetic distributions of thermal f,{R,y) and non- 
thermal f n t(R, y) electron components for distances of R — 50r g 
(first panel) and R = \0r g (second panel), respectively, for 
A = 10. 



Taking into account the fact that adiabatic compression in- 
creases the perpendicular electron momentum only (and so we 
may take sin iff « 1 for the pitch-angle) and convolving with the 
electron distribution J34> . the final expression for the accretion 
flow emission spectrum is 




10'* 10 la 10'° 10" 10'° 10' 3 10 £u 10" 10" 

Frequency 

Figure 5. Decomposition of a single black hole (with the mass 
1OM ) emission spectrum into thermal and nonthermal parts. 
The accretion rate is 1.4 • 10 10 g/s, which corresponds to m — 
1(T 8 , and £ = 0.1. 




Frequency 



Figure 6. Accreting 10M o black hole spectra for accretion rates 
from m = 10~ 9 till m = 10~ 5 (M = 1.4 ■ 10 9 g/s and M = 
1.4 • 10° g/s, correspondingly), for£ = 0.1. 



5. Luminosity 

The luminosity of the accretion flow may be computed by inte- 
grating expression d63i (Shapiro 1973a, Ipser & Price 1977 1 

CO COS &* CO 



L v cc J R- 5 ' A dR ff(R,y) jf YZ-^^F^y cos 6dy.(63) L = 8 ^ J RldR j j Jy' dv ' 

1 1 -1 1-10 

x(i-— ) (i -yS 2 ) 3/2 (i - p cos ey 3 d cos e. (64) 



Fig- |U shows a decomposition of the accretion flow spec- 
trum into thermal and nonthermal parts. The fraction of the 
nonthermal emission is small in the optical range, but domi- 
nates in harder spectral bands. The spectrum becomes flatter 
with the accretion rate decrease due to the increase of the elec- 
tron radiative loss timescale. The high-energy spectral cut-off 
is determined by radiative energy losses and by the upper limits 
of the accelerated electron gamma-factor (see Fig|6jl. 



Dividing by Mc 2 and taking into account the expression for 
synchrotron luminosity of a single electron j38> . we get the ex- 
pression for the efficiency of the thermal component emission 



It = 4m j ^K(R)dR, 
1 



(65) 
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where 



K(R) 



1 f(l - l/R) l/2 (l -p 2 f 12 



(1 -pxf 



dx. 



(66) 



At low accretion rates (m «: 10 5 ) a reasonably good ap- 
proximation for this quantity is 



rj, = 4.5 • 10 4 m, 

which corresponds to the thermal luminosity of 
L, = 9.6 ■ 10 33 M 3 l0 n 2 (V 2 + c])\l erg/s. 



(67) 



(68) 



The efficiency of the nonthermal component emission may 
be easily estimated from the radial distribution of the dissi- 
pated magnetic energy dE/dRdt, its fraction £ carried out by 
accelerated particles and the evolution of its mean square of 
the gamma-factor during the fall towards the horizon y\{R) as 

( fFvi \ 

K(R)dR, (69) 



" 4 J r 2 (R )R 2 J { dR I 



r 2 (R )R 2 

1 1 

where the mean radiative energy losses of a single electron are 
~7 



(70) 



The evolution of the mean square of gamma-factor y 2 b (R) of the 
beam is given by expression (IA.4t derived in AppendixlAl 

So, for the efficiency of the non-thermal component emis- 
sion we have 



Vnt = 



~\2 



CdR f dR 



K(R) 



xS 2 



AT(R ) 



R 



1 - 



R 



13/8\ 



(71) 



The results of numerical computations according to these 
formulae may be seen in Fig0 

6. Properties of flares 

We now discuss some temporal properties of the accretion flow 
emission. Neglecting the complex spatial structure of the ac- 
cretion flow we may consider the thermal electron component 
emission Lo as a constant background with a highly variable 
nonthermal emission component superimposed. The variabil- 
ity of the latter is mostly due to the discrete nature of parti- 
cle acceleration in the current sheets. Assuming for simplic- 
ity that each non-thermal electron accelerates only once and 
then free-falls and evolves due to adiabatic heating and radia- 
tive losses only, we get for the light curve of a single flare (i.e. 
a beam of N electrons ejected by a reconnection event at some 
distance Rq with distribution function fl23i and mean gamma- 
factor r(Ro) = At{Rq)) an expression 



AL = ^o- T c (y 2 b - 12t 2 ) ^NK(R), 



(72) 




Accretion rate 



Figure 7. Efficiencies of the synchrotron emission of thermal 
and non-thermal electron components of the accretion flow. 

where the quantity y 2 b is derived in Appendix |A| and repre- 
sents the mean square of the gamma-factor of the electron 
beam ejected at Ro and free-fallen to R. The mean square of 
the Maxwellian distribution d36i gamma-factor value of 12r 2 
is also taken into account (we need to subtract this term since 
each electron must occur in the total luminosity expression only 
once). The K(R) coefficient describes relativistic effects of the 
radiation reduction. 

By introducing an effective volume of the current sheet 
V = N/n(Ro) and dividing by the mean luminosity of the non- 
thermal component L, = r] n ,Mc 2 , we get 



AL 



llnrjntR^ 2 



(Y\ y 2 b -\2T 2 

R5/2 



K(R). 



(73) 



Transition to the light curve may be performed by substituting 
the dependence of distance on time for the falling gas element 

as 

\2/3 



3a c 



(74) 



The temporal structure of the flare is complex (see HiglHJ. It 
consists of several parts with the domination of different phys- 
ical processes on the electron cloud evolution. The first stage 
is the acceleration of the electrons in the current sheet (see 
section l2~4l . which occurs on the Alfvenic time scale, but the 
exact temporal structure depends on the properties of the cur- 
rent sheet and thus is beyond the scope of the article. The sec- 
ond stage is the fast emission decay due to synchrotron energy 
losses on the corresponding energy scale; it lasts until the es- 
tablishment of the equilibrium between energy losses and adi- 
abatic heating rates. The third stage is generally the longest 
one, at least for Rq < 10 3 , and has the time scale of a free-fall. 
Finally, near the event horizon the relativistic effects (time con- 
traction, gravitational redshift, Doppler effect and the photons 
capture by the BH) prevail and flare emission decays on the 
r g /c time scale as the electron cloud approaches the horizon. 
This stage is very important as it directly reflects the metrics 
near the black hole. 
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Table 1. Thermal and nonthermal luminosity, magnitudes, fluxes and flares rate of a 10 M black hole in a various ISM conditions 
at a 100 pc distance. 



n V m L, L„, m v m% EPIC 1 HRC 2 JEM-X 3 IBIS 4 XRT 5 Flare rate' 

[cm -3 ] [km/s] [erg/s] [erg/s] [erg/s/cm 2 ] [erg/s/cm 2 ] [erg/s/cm 2 ] [erg/s/cm 2 ] [erg/s/cm 2 ] [10 3 /s] 
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1 XMM EPIC band flux. Sensitivity is 10~' 4 erg/s-cm 2 for typical survey exposures. 

2 Chandra HRC band flux. Sensitivity ranges from 10~ 14 erg/s-cm 2 to 10~ 16 erg/s-cm 2 (deep surveys) 

3 INTEGRAL lem-X band flux. Sensitivity is 4 • 10~" erg/s-cm 2 for 100 ksec exposure 
3 INTEGRAL IBIS band flux. Sensitivity is 10~' 5 erg/s-cm 2 for 1000 ksec exposure 

5 Swift XRT band flux. Sensitivity is 5 ■ 10~ 13 erg/s-cm 2 for 1 ksec exposure 

6 Flare rate is for the flares with mean luminosity greater than 0. 1 % of mean nonthermal one 



ejection multiplied by the speed of light) is smaller than that of 
the magnetic field loops, then we see the beam emission only 
when it is pointed towards us. The timescale of such flares is 
much shorter than the free-fall one, and they can be detected 
only in high time resolution observations. On the other hand, 
their properties reflect the magnetic field structure, and so its 
search is very important. 

A significant amount of nonthermal synchrotron emission 
also falls in the optical range (see Table ^ an d Fig ll li . This 
means that short flares like that shown in Fig|9]may be detected 
in observations with a 1 fis time resolution at the 6 meter tele- 
scope in the framework of the MANIA experiment (Beskin et 
al. ri997l . 



7. Discussion 

The analysis of existing data on possible black hole masses 
and velocities is performed in comparison with the interstel- 
lar medium structure. It is shown that in the majority of cases 
in the Galaxy (> 90%), the accretion rate m = Mc 2 jL e dd can- 
not exceed 10~ 6 - 10~ 7 (see FigQ. The basis of our analysis 
is the assumption of energy equipartition in the accretion flow 
of Shvartsman dl971> . Currently, there have been no successful 
attempts to build a magnetized plasma accretion theory without 
this assumption - the work of Kowalenko & Melia ( 1 999 1 led to 
unphysical results, while Scharlemann ( 1983) has some math- 
ematical errors in computing the back-reaction of the magnetic 
field on the matter. 



Fig|9] shows the sample light curves of flares with Rq = 5, 
R = 10, R = 20 and R Q = 30, for the accretion rate m = 10" 8 
and the equivalent volume V — r*. The unknown first stage 
(flare front) is excluded. 

By averaging over the flare light curve we may get the de- 
pendence of mean flare amplitude < AL/L„, > on the ejection 
radius Ry. This is shown in Fig^| It is clear that the strongest 
flares are produced by reconnections very near the event hori- 
zon (at several r g ). The rate of a flares may be easily computed 
as 

dX dE £ 



?o) J \dR L 



K(R)dR (75) 



dR dt dR dt ET(R 
where E is the total energy of a flare as is seen at infinity 
E = L„, j < AL/L nt > dt (76) 

By integrating this expression over the range of Rq where 
< AL/L„ t > is greater than some threshold value we may get 
the rate of flares with given amplitudes. Results of such compu- 
tation for the rate of flares with an amplitude greater than 0.1% 
of the mean nonthermal luminosity (< AL/L nt > > 10~ 3 ) are 
given in the Tabled 

The fine temporal structure of these flares may be related 
to the motion of the electron beam in the magnetic field (this 
motion, as it has been noted before, is finite; electrons do not 
leave the small volume with the size determined by the Larmor 
radius and magnetic field line topology). So, if the character- 
istic size of the beam (the characteristic timescale of particle 
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Adiabatic heating 



Free-fall towards the horizon 
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Figure 8. Internal structure of a flare as a reflection of the 
electron cloud evolution. The prevailing physical mechanisms 
defining the observed emission are denoted and typical dura- 
tions of the stages are shown. 
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Figure 10. Mean amplitude of a flare as a function of ejection 
radius Rq. 




time until event horizon 



Figure 9. Light curves of separate non-thermal flares - beams 
of accelerated electrons ejected at Ro - for m = 1(T 8 . Time is 
measured in units of r g /c, luminosity in units of total thermal 
luminosity. The stage of particle acceleration (which defines 
the flare front) is omitted as its shape depends on an unknown 
temporal structure of the particle acceleration process, which is 
beyond the scope of our consideration. 



The accreting plasma is initially collisionless, and it re- 
mains so until the event horizon. The electron-electron and 
electron-ion free path A ~ 2.4 ■ 10 3 r 2 n~' even at the capture 
radius is as high as ~ 10 12 cm. Only the magnetic fields trapped 
in plasma (the proton Larmor radius at r g is 10 cm) make it pos- 
sible to consider the problem as a quasi-hydrodynamical one; it 
is only due to the magnetic field that the particle's momentum 
is not conserved, allowing particles to fall towards the black 
hole. In addition, the magnetic field effectively "traps" parti- 
cles in a "box" of variable size, which allows us to consider its 
adiabatic heating during the fall; a correct treatment of such a 
process shows that for magnetized plasma such heating is 25% 
more effective than for ideal gas (see Sec B.U . Therefore, the 
plasma temperature in the accretion flow grows much faster and 
electrons become relativistic earlier - R re i « 6000 in contrast 
to R re { » 1300 in Bisnovatyi-Kogan & Ruzmaikin ( 1974 1 and 




Figure 11. Contribution of nonthermal component to the opti- 
cal emission /L° pt as a function of accretion rate 



R re i ~ 200 in Ipser & Price ( 1982 1. The accretion flow is much 
hotter, and our estimation of "thermal" luminosity 

L = 9.6 • 10 33 M 3 « 2 (V 2 + c 2 )r 3 erg/s (77) 

is significantly higher than those of Ipser & Price ( 119821 

Ljp = 1.6 • 10 32 M 3 » 2 (V 2 + c 2 )[ 3 erg/s (78) 



and Bisnovatyi-Kogan & Ruzmaikin ( 1974 
Lbkr = 2 • 10 33 M 3 « 2 (V 2 + c 2 ): 3 erg/s, 



(79) 



while the optical spectral shape is nearly the same. 

We considered dissipation of the magnetic energy in the 
turbulent current sheets (Pustilnik 1997) as a mechanism sup- 
porting equipartition. In this process, the electrons are ejected 
with plasma from the current sheet and are accelerated. These 
electrons have a power-law energy distribution and its emission 
spectrum is flat up to the gamma band (Figs|5]and|6l>. An im- 
portant property of the nonthermal emission is its flaring nature 
- the electron ejection process is discrete; typical light curves 
of single beams are shown in Fig.|S]and|5] 
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Tablensummarizes the observable parameters of accreting 
black holes. 

The black holes have significant luminosity only when 
they are located in the warm hydrogen regions that occupy 
about 50% of the Galaxy volume (McKee & Ostriker[T977}. 
In such cases the dimensionless accretion rate lies in the 1CT 8 
- 10~ 6 range, and the bolometric luminosity of the accretion 
flow is 3 ■ 10 28 -10 33 erg/s, depending on the BH velocity and 
the fraction of magnetic energy transformed into the motion 
of electrons (see Table As a result, the black hole at a 
100 pc distance (a sphere with this radius must contain sev- 
eral tens of such objects, see Beskin et al ( 2000 1 and Agol & 
Kamionkowski (2002 1) looks like a 15-25™ optical object (due 
to the "thermal" spectral component) with a strongly variable 
companion in high-energy spectral bands ("nonthermal" com- 
ponent). The hard emission consists of flares, the majority of 
which are generated inside a 5r g distance from the BH (see 
Figs. l8l9fL0l . These events have the durations ~ r g /c (~ 10~ 4 
s), a rate of 10 3 -10 4 flares per second, and an amplitude of 2%- 
6%. As it is seen from Tabled the BH variable X-ray emission 
may be detected by modern space-borne telescopes. 

Optical emission consists of both a quasistationary "ther- 
mal" part and a low-frequency tail of nonthermal flaring emis- 
sion. The rate and duration of optical flares are the same as 
X-ray ones, while their amplitudes are significantly smaller. 
Indeed, the contribution of nonthermal component to the op- 
tical emission (see Fiell 1> is approximately 2 • 10~ 2 for m = 
10~ 8 - 10~ 6 , so the mean amplitudes of optical flares are 0.04%- 
0.12%, while the peak ones may be 1.5-2 times higher and 
reach 0.2%. Certainly, it is nearly impossible to detect such sin- 
gle flares, but their collective power reaches 18-24™ (see Table 
Q and thus may be detected in observations with high time res- 
olution (< 10~ 4 s) by the large optical telescopes. 

Our separation of the accretion flow emission into station- 
ary "thermal" and flaring "nonthermal" components is a rough 
approximation made for the estimation of the qualitative pic- 
ture of BH observational appearances. The real behaviour of 
the accreting plasma may be much more complicated. 

The accretion flow is a complex dynamical system with 
nonlinear feedback. This is ensured by the plasma oscillations 
generated in each reconnection event, with beams of acceler- 
ated electrons and clouds of magnetized plasma ejected from 
current sheets. All these agents may act as triggers for already 
"prepared" inhomogeneities which turn on magnetic energy 
dissipation processes. This situation seems to be similar to the 
Solar one which determines its flaring activity, and also to the 
case of UV Cet stars and maybe accretion disks of X-ray bi- 
naries and active galactic nuclei. All these non-stationary pro- 
cesses are characterized by power-law scalings of flare energies 
with similar slopes of 1 .5-2 at a very wide range of energies - 
from 10 23 erg/s for the Sun to 10 45 erg/s for quasars. The uni- 
versality of these processes may be interpreted in the frame- 
work of a fractal approach as done by Bak, Tang & Weisenfeld 
(TTW71. Lu & Hamilton (1TWT1 . Lu et al. (TH?33l . Anastasiadis, 
Vlahos & Georgoulis l l 19971 . Kawaguchi & Mineshige ( 1999 1, 
Pustilnik d!997i . This means the realization (at least in ac- 
tive phases) of some collective state, sometimes called "self- 
organized criticality" (Bak, Tang & Weisenfeld 1987 1), which 




1 10 100 

R 

Figure 12. Fraction of thermal synchrotron emission that 
comes from inside a given radius R 

is characterized by the same behaviour of the parameters on all 
scales. These are percolation processes. There is evidence that 
accretion flow is in this state, and so its observational appear- 
ance (at least that related to non-stationary processes) may be 
predicted and interpreted in the framework of this approach. 
Initial steps in this direction have been made by Beskin & 
Karpov ( 2002 ), but need to be refined. 

8. Conclusions 

During recent years the number of works dealing with single 
stellar-mass black holes has significantly increased. Some are 
purely theoretical (Punslv fT9~98al 1 1 998b| Gruzinov & Quataert 
1999, Abramowicz et al. 12002 1 and other provide discussions 
of their observational detection (Heckler & Colb 1996; Fujita 
et al. [T99l Beskin et al. l2000l Agol & Kamionkowski I2UU1 
Chisholm et al. 2003 1. The importance of experiments in strong 
gravitational fields has been recently noted by Damour ( 1998 1. 
Kramer et al. (2004i discussed the new possibilities of black 
holes metric study by the investigation of radio pulsar - BH 
binary systems with new generation of radio telescopes. 

In this work we tried to concretize physical properties 
of plasma accreted onto the black hole within the classical 
paradigm of equipartition of Shvartsman J197U . Assuming the 
discrete nature of the magnetic energy dissipation processes in 
current sheets allows us to clarify the shape of the synchrotron 
spectrum of the accretion flow. A hard highly non-stationary 
nonthermal spectral component appears as an emission of ac- 
celerated particles. The beams accelerated in the current sheets 
can generate very short flares, providing information about the 
neighborhood of the event horizon (Fig|9). On the other hand 
it is clear from Fig [R] that at low accretion rates a significant 
amount of thermal synchrotron radiation is generated inside 3r g 
- this means that the behaviour of this component will reflect 
the properties of space-time in strong gravitational fields too. 

It is clear that the search for a black hole strategy may be 
modified in accordance with such results. Optical high time 
resolution studies of X-ray sources may be very important. 
Single black holes may be contained inside the known sta- 
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tionary gamma sources (Gehrels et al. 2000) as well as objects 
causing long microlensing events (Paczynski 1991 1. Thus it is 
very important to look for X-ray emission as well as for fast 
optical variability of these objects. Sample observations of the 
longest microlensing event MACHO 1999-BLG-22 (Bennett et 
al. 2002 1, a stellar-mass black hole candidate, have been per- 
formed at the Special Astrophysical Observatory of RAS in the 
framework of the MANIA experiment in 2003-2004 (Beskin et 
al. l2005l . 

The best evidence will be provided by the synchronous high 
time resolution observations in optical and X-ray ranges. 

Detection of the event horizon signatures cannot result from 
statistical studies. A detailed study of each object is needed to 
detect its specific appearance. 
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Appendix A: Moments of single accelerated beam 
energy distribution 

For the case of the electron beam ejected from the current sheet 
at some radius Rq and free-falling to R not interacting with the 
background particles and beams ejected below the moments of 
the energy distribution may be written using (I46> as 



n(R) = 



n(R) = 



I CO 

J f(Ro, 
*- 1 



yo)dyo 



J f(R ,y )dy J 

VI / 1 



f f( R o> yo)yodyo 

J C l (R,R )yo + C 2 (R,R ) 
l 

f(Ro, yo)y a dyi) 



CdR,Ro)yo + C 2 (R,Ro) 



(A.2) 



For the initial distribution of the form A23I this may be eval- 
uated as 



y b (R) = T(R m J — 



y 2 b (R) = r 2 CRo)S2 



where new functions 




3/2\ 



3/2\ 



MA,B) = 

= i(A-j?.« ap (|)m 1 (|)) 
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(A.v+B) 2 



.(-A + (1+B).exp(f ^(f)) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



are introduced and an expression for the integral exponent 

CO 



Ei„(x) 



is used. 



J ~F 



(A.7) 
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